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Overview 


AppUcation  of  CFD  to  aircraft  stability  &  control  (S&C) 

•  Background 

•  Air  Force  SEEK  EAGLE  Office  (AFSEO) 

•  The  problem... 

•  Suggested  approach 

•  Examples: 

■  static:  F-16XL  (CAWAPI)  /  F-18 

■  Maneuvering:  F-16C 

•  Conclusion 
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Background 


•  Virtually  every  new  a/c  and  armament 
program  encounters  unexpected  S&C 
problems 

■  Current  S&C  modeling  techniques  inadequate 

■  Complex  S&C  issues  require  nonlinear,  CFD- 
in-the-loop  simulations  with  control  surfaces 
modeled 

•  New  CFD  methods  combined  with  HPC 
resources  ... 

■  address  S&C  problems  that  previously 
were  untenable 

■  overcome  experimental  gaps 

■  complement  experiments 

■  result  in  safer  flight  testing 

■  reduce  amount  of  wind  tunnel  and/or 
flight  testing 
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Computational  S&C  Approaches 


•  "Brute-force"  approach  to  filling  an  aerodynamic 

database  (Murman  eta/,  Rogers  et  al) 

■  30,000  solutions  256  procs/16  mill  cells  158  years 

■  Only  static  data  at  discrete  points 

•  Data  fusion  approach  (Tang  eta/) 

■  Intelligent  agents  to  combine  low/high  order  solutions 

•  Automated  Cartesian-based  methods  (m  urman) 

•  Reduced-frequency  approach  to  dynamic  derivatives 

(Murman) 

•  Many,  many  others... 

•  New  RTO  task  group  ^  AVT  166 
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USAF  SEEK  EAGLE  Office  (AFSEO) 


•  The  SEEK  EAGLE  program  is  the  standard  for  the  aircraft-stores 
certification  process  for  the  US  Air  Force 

•  Provides  Quick  Reaction  Certifications  (QRC),  Certification 
Recommendations  (CR),  and  Flight  Ciearances  (FC) 
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The  AFSEO  Problem 


1997  2002  2007 
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...and  Other  Issues 


•  F-16  and  F-22  ATLAS  aerodynamic  database  limitations 

■  Limited  in  the  number  of  configurations 

■  Only  symmetric  configurations  available 

■  No  data  for  modern  stores  and  suspension  equipment 

•  Unable  to  predict  where  instabilities  occur  in  the  flight 
envelope  (if  at  all) 

•  Long  history  of  flight  test  programs  where  lots  of  $$ 

spent  to  find . absolutely  nothing 
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Suggested  Approach 


•  Focus  on  specific  flight  conditions  (if  possible) 

•  Simulate  closed-loop,  full-scale  a/c  at  critical  conditions 
with  a  single,  complex  and  efficient  maneuver  (possibly 
non-flyable) 

•  Generate  nonlinear,  dynamic  reduced-order  models  for 
aerodynamic  loads 

•  Use  model  for  S&C  analysis,  flight  simulation,  control 
system  design,  etc. 

■  Continuous  data 

■  Derivatives  computed  analytically 

■  Allows  SEEK  EAGLE  Office  flexibility  to  handle  any  new 
configuration  and  independence  from  contractors 
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Suggested  Approach 


1)  Model  Training 
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3)  Model  Prediction 
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Examples 


Flow  Solver  -  Cobalt 


Numerical  Modeling 

•  Unstructured,  cell-centered, 
finite  volume  CFD  code 

•  Hybrid  grids  (arbitrary 
elements),  overset  grids 

•  Space  discretization: 

■  Godunov's  first-order  accurate 
exact  Riemann  solver 

■  Second-order  accuracy  through 
least-squares  reconstruction 

•  Time  discretization: 

■  Point-implicit  with  Newton  sub¬ 
iterations  for  time  accuracy 

•  6DOF  and  ALE  for  rigid  body 
motion 


Physical  Modeling 

•  Compressible  Euler  and  NS, 
laminar  and  RANS 

•  Equilibrium  air  physics 

•  Turbulence  models: 

■  Spalart-Allmaras  (SA) 

■  Mentor's  SST 

■  Wilcox's  1998  k-(o 

■  Detached  Eddy  Simulation  (DES) 
for  S-A  and  SST 

■  Curvature  corrections  (SARC) 
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F-16XL 


•  RTO  AVT-llS  (CAWAPI) 

•  Comprehensive  validation  of  CFD  methods  against 
available  flight  test  database 


•  Time-accurate  DES 
at  full-scale  flight 
Reynolds  numbers 


calculations 


Morton,  S.A.,  McDaniel,  D.R.,  and  Cummings,  R.M.,  ''F-16XL  Unsteady  Simulations  for 
the  CAWAPI  Facet  of  RTO  Task  Group  AVT-llS,"  AIAA  Paper  2007-0493,  Jan.  2007. 
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F-16XL  (CAWAPI) 


F-16XL  SARCDES  Solution  at  FC7 

Alpha  =  11.89  deg  Mach  =  0.304  Re#  =  44.4E8 

Vorticity  Magnitude  Iso-Surface  (250  1/sec)  Colored  By  Pressure  (psi) 
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F-16XL  (CAWAPI) 
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F-16XL  (CAWAPI) 
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F-18C  (HARV) 


•  High  Angle  of  Attack  Research  Vehicle  (HARV)  -  NASA-Dryden 

•  Rich  source  of  flow  viz,  surface  pressures,  and  aeroelastic  data 

•  Numerical  investigation  to  try  and  accurately  predict  flow 
impinging  on  vertical  tails  (M=0.28,  Re#  =  13E6,  a=30  deg) 


F-18  HARV 

Smoke  Test 
late  1980's 

Dry  den 

Flight  Research  Center 


Morton,  S.A.,  Cummings,  R.M.,  and  Kholodar,  D.B.,  "High  Resolution  Turbulence 
Treatment  of  F/A-18  Tail  Buffet,"  AIAA  Paper  2004-1676,  April  2004. _ 
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Play 
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Power 


F-18C  (HARV) 


Breakdown  Pressure 


Strouhal  number,  fc/U 

^  00 


22  June  2007 


17 


x^l  Breakdown 


F-18C  (HARV) 
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Full-Scale  F-16C  Motion  Simulations 


Half-span  grid  with  3.4  x  10® 
cells  and  prismatic  layers 

Cells  concentrated  in  the 
St  rake  vortex 

Forebody  bump,  diverter, 
ventral  fin  modeled 

Engine  mass  flow  modeled 

Flow  conditions: 

■  M^=  0.25/0.60 

■  Re=14.7  X  10V42  x  lO® 

Numerical  parameters: 
■At=0.0004s/0. 0002s  (At*=0.01) 

■  5  Newton  sub-iterations 

■  DES  based  on  SA  with  RC 
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a-Sweeps  &  Dynamic  Scale  Effects 


0  5  10  15  20  25 

a  (deg) 


•  Angle  of  attack  sweeps  at  various  pitch  rates 

-  60  deg^sec  (1^3  scale,  dr=0.01 ,  Re/f^2E6) 

120  deg/sec  (1/3  scale,  dr=o,oi,  RemF2E6J 
-  30  deg^sec  (full  scale,  dr^0.03,  R©'ft=  3.7ES)) 
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Aircraft  System  Identification  (SID) 


•  Normally  applied  to  flight-test  data  to  generate 
aerodynamic  models  (assumed  time-invariant) 

•  Can  use  SID  techniques  to  analyze  CFD  data 
computed  for  aircraft  in  prescribed  motion 

•  System  Identification  Programs  for  Aircraft 
(SIDPAC)  by  Klein/Morelli  at  NASA-Langley 
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Pitch  Chirp  Training  Maneuver 
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Pitch  Chirp  (X=1.0) 
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SID  of  Pitch  Chirp 


1.  Least-squares  linear  model:  =Cq  +  q  +  q 


2.  Nonlinear  multivariate  polynomial  model :Cjy  =Cg  +Cja  +  C2q  +  Cjq  +  ... 

•  Model  structure  determined  using  +c,a^  +C,aq^  +Cjaq  +  ... 

multivariate  orthogonal  functions 

generated  by  Gram-Schmidt  +Cgqq  +  Cgq^q  +  C^Qq^  +Cnq^ 

orthogonalization,  ordered  by  dynamic 
programming  „,2 

•  Retained  modeling  functions  expanded 
into  ordinary  multivariate  polynomial  <" 

o  0  09 

•  Parameters  estimated  using  maximum  - 

likelihood  technique  f” 

^  0.07 
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System  ID  of  Pitch  Chirp 
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static  Validation  vs.  CFD/ATLAS 
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Dynamic  Validation  vs.  CFD 


model  prediction  for  validation  data  (nord=4,  maxord=4) 


validation  data  (CFD) 

multivariate  polynomial  model,  linear  model  fit  to  estimate  sig2 
multivariate  polynomial  model,  constant  model  fit  to  estimate  sigS: 
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Dynamic  Validation  vs.  ATLAS 
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Dynamic  Validation  vs.  ATLAS 
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1  Hz  pitch  osciiiation  about  a  =  15  deg 
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Angle  of  Attack  (deg] 


Yaw  (X=1.5) 


Play 
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Composite  Pitch-Yaw  Chirp 


•  Single  motion  input  to  create  a  model  including  motion  about  two 
axes  ^  a  =  15+15  deg,  p  =  0+15  deg 


•  Input  signals  are  made  orthogonal  by  setting  Xto  1.0  for  pitch 
and  varying  X  until  dot  product  of  the  two  signals  is  zero  resulting 
in  X  of  1.47  for  yaw  signal 


Qr  +  Cyf^aP^  +  Cj,«^  +  P  +  q  +  +  C\(,P^ 
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Composite  Pitch-Yaw  Chirp 
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Rotation  About  Axis  (deg) 


Pitch  Doublet  Flight  Test  Maneuver 


•  Gathered  actual  strip  chart  data  from  a  flight  test 

•  Created  the  motion  file  that  forces  the  F-16C  through 
translations  and  rotations  giving  similar/same  behavior 


•  Predict  maneuver  with  reduced-order  model  and 

compare  WbclelPrecicBonlbrCL 
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Pitch  Doublet  Flight  Test  Maneuver 

Simulation  run  at  M=0.6  and  h=5,000  ft  with  full  span  grid 
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2.5g  Wind-up  Turn  Flight  Test  Maneuver 

Simulation  run  at  M=0.6  and  h=5,000  ft  with  full  span  grid 
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Conclusions  and  Ongoing  Work 


state-of-the-art  CFD  techniques  capable  of  predicting  static,  unsteady 
flow  fields  with  complex  aerodynamic  behavior 

Motion+modeling  technique  shows  promise  and  has  many  benefits 

Goal  is  to  direct/complement  experimental  techniques  -  not  replace  them 

Need  good  full-aircraft  S&C  experimental  data  for  validation  (AVT-166) 


Ongoing  Work; 


More  accurate  F-16C  grid  with  stores  and  pylons 

Comparison  of  CFD  and  model  predictions  with 
ATLAS/flight  test  data 

Integration  of  moving  control  surfaces  and 
inclusion  into  the  SIDPAC  modeling  process 

Investigation  of  range  of  validity  of  models 
across  flight  envelope 


F-16C  Trailing  Edge  Flaperon,  ±30  deg 


Aeroelastic  deformation 
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